home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / modelers / geomview / source.lha / Geomview / src / lib / gprim / bezier / bezbound.c next >
Encoding:
C/C++ Source or Header  |  1992-10-01  |  1.7 KB  |  62 lines
  1. /* Copyright (c) 1992 The Geometry Center; University of Minnesota
  2.    1300 South Second Street;  Minneapolis, MN  55454, USA;
  3.    
  4. This file is part of geomview/OOGL. geomview/OOGL is free software;
  5. you can redistribute it and/or modify it only under the terms given in
  6. the file COPYING, which you should have received along with this file.
  7. This and other related software may be obtained via anonymous ftp from
  8. geom.umn.edu; email: software@geom.umn.edu. */
  9.  
  10. /* Authors: Charlie Gunn, Stuart Levy, Tamara Munzner, Mark Phillips */
  11.  
  12. #include "bezierP.h"
  13.  
  14. BBox *
  15. BezierBound( bezier, T )
  16.     register Bezier *bezier;
  17.     Transform T;
  18. {
  19.     register float *p;
  20.     register int v, n;
  21.     HPoint3 min, max;
  22.  
  23.     /* only support 3 and 4 dimn patches */
  24.     if(bezier->dimn == 4) {
  25.     /* It may not make sense to find the bbox of these 4-D control points.
  26.      * Dice it into a real mesh and take the bbox of that.
  27.      */
  28.     if(bezier->flag & BEZ_REMESH ||
  29.        bezier->mesh == NULL || bezier->mesh->p == NULL) {
  30.         if(BezierReDice(bezier) == NULL)
  31.         return NULL;        /* Oh no */
  32.     }
  33.     return MeshBound( bezier->mesh, T );
  34.     }
  35.  
  36.     if (bezier->dimn != 3) {
  37.        GeomError(0,"BezierBound: invalid dimension %d",bezier->dimn);
  38.        return(NULL);
  39.     }
  40.  
  41.     n = (bezier->degree_u + 1) * (bezier->degree_v + 1);
  42.  
  43.     p = bezier->CtrlPnts;
  44.     Pt3Transform( T, p, &min );
  45.     min.w = 1;
  46.     max = min;
  47.     for( v = n; --v > 0; ) {
  48.     Point3 t;
  49.  
  50.     p += 3;
  51.     Pt3Transform( T, p, &t );
  52.     if(min.x > t.x) min.x = t.x;
  53.     else if(max.x < t.x) max.x = t.x;
  54.     if(min.y > t.y) min.y = t.y;
  55.     else if(max.y < t.y) max.y = t.y;
  56.     if(min.z > t.z) min.z = t.z;
  57.     else if(max.z < t.z) max.z = t.z;
  58.     }
  59.  
  60.     return (BBox *) GeomCCreate (NULL, BBoxMethods(), CR_MIN, &min, CR_MAX, &max, NULL);
  61. }
  62.